Strong pressure dependence of the magnetic penetration depth in single crystals of 
the heavy fermion superconductor CeCoIns studied by muon spin rotation 
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In the tetragonal heavy fermion system CeCoIns the unconventional superconducting state is 
probed by means of muon spin rotation. The pressure dependence (0—1 GPa) of the basal-plane 
magnetic penetration depth (\ a ), the penetration depth anisotropy (7 = A c /A a ) and the temperature 
dependence of 1/X 2 (i = a, c) were studied in single crystals. A strong decrease of A a with pressure 
was observed, while 7 and A?(0)/A?(T) are pressure independent. A linear relationship between 
1/Aa(270 mK) and T c was also found. The large decrease of A a with pressure is the signature of an 
increase of the number of superconducting quasiparticles by a factor of about 2. 



Unconventional superconductors are characterized by 
their proximity to different instabilities. In heavy fermion 
systems superconductivity is often found in the region of 
the phase diagram where a weak magnetic phase disap- 
pears [H. However, in some systems superconduc- 
tivity is also detected in proximity of a valence phase 
transition. 

CeCoIns is a prototypal heavy fermion superconduc- 
tor 0] at the focus of numerous studies owing to the 
proximity of quantum criticality. This proximity is re- 
flected by the pronounced non-Fermi liquid features Q 
and the highest superconducting (SC) transition temper- 
ature T c — 2.3 K [3| among the Ce based heavy fermions. 
In addition, this tetragonal system is characterized by a 
uasi two-dimensional Fermi surface [|1 Q and a two-gap 
unconventional SC state with d-wave symmetry [8|. 
In order to clarify the relation between the SC phase 
and quantum criticality, the evolution of basic SC pa- 
rameters with respect to a tuning parameter is required. 
In this letter the pressure (p) and temperature (T) de- 
pendence of a fundamental SC quantity — the magnetic 
penetration depth (A*) — was studied. Here i — a,c cor- 
responds to a screening current flowing along the main 
crystallographic directions: perpendicular, respectively 
parallel to the c-axis. A^ is obtained through the precise 
magnetic field distribution in the SC vortex state probed 
by transverse-field (TF) muon-spin rotation (/iSR). 

CeCoIns is in the clean limit with a ratio of coherence 
length to mean free path £/l < 0.02 (£ < 8.2 nm) 
In this limit the penetration depth can be written in the 
London model as: 



lAi = ti e 2 n s /m* 



(1) 



Here /zo is the vacuum permeability, e the electron 
charge, n$ the number density (number of supercon- 



ducting quasiparticles), and m* the effective quasipar- 
ticlc mass. 

The /iSR experiments were performed at the Swiss 
Muon Source (S/iS), Paul Scherrer Institute (PSI), 
Switzerland, using the GPD (under p) and LTF (ambient 
p, low T) spectrometers. In a TF-/xSR experiment spin 
polarized positive muons are implanted into a sample in 
an external magnetic field [IqH (field cooled from above 
T c ) applied perpendicular to the initial muon-spin polar- 
ization. In the presence of a magnetic field at the muon 
site the muon spin precesses at its Larmor frequency 
^ti = li+B^ (7^ = 8.516 • 10 8 rad s _1 T _1 is the gyromag- 
netic ratio of the muon) before decaying with a life time 
of — 2.2 fis into a positron and two neutrinos. Due 
to parity violation the decay positron is preferentially 
emitted along the muon spin direction. Forward and 
backward positron detectors with respect to the initial 
muon polarization are used to monitor the /xSR asym- 
metry spectrum A(t). 

Single crystals of CeCoIn5 were grown by indium flux 
method (rare earth from [l(j), centrifuged and etched 
in HC1 solution to remove the indium excess. Thin plate- 
like single crystals were obtained with their large faces 
corresponding to the (001) basal plane. Using this par- 
ticular geometry, two samples were prepared, consisting 
of ~10, respectively ^200 crystals glued together with 
G.E. varnish, as sketched in Fig. [21 The mosaic sample 
(c-axis normal to the plane) was studied with the LTF 
spectrometer. The cylindrical-like sample (a-axis is the 
main axis) was mounted in a piston cylinder pressure cell 
of CuBe alloy with Daphne oil as a pressure transmission 
medium [ll[ and measured with the GPD spectrometer. 
The actual pressure in the cell was determined by the T c 
of a small piece of indium. 

For different pressures an angular scan consisting of 
5-8 TF-/1SR spectra was taken at T ~ 270 mK with an 
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FIG. 1. (Color online) /J.SK asymmetry spectrum A(t) of 
CeCoIns at ambient pressure in the normal (black squares) 
and SC state (red circles) taken with the LTF spectrometer 
at fioH = 50 mT for H \\ c. The data are shown in a frame 
rotating at a frequency corresponding to 45 mT p3 |. Solid 
lines represent fits to the data (see text). 



applied field hqH ~ 50 mT forming an angle 9 with the 
sample's crystallographic c-axis. For p = GPa, 0.2 GPa, 
and 0.6 GPa a temperature scan was also recorded for 
9 = 0° (H || c) and 9 = 90° (H _L c). The field was 
chosen to be higher than the critical field of bulk indium 
{fJ.oH C 2(0) = 23 mT) to avoid artifacts due to possible 
residual flux from the growth. For comparison the val- 
ues of the Pauli limited critical fields for CeCoIns are: 
H H cl ~ 10 mT (H || c and H _L c) Q, n H c2 = 5 T 
(H || c) and 11.5 T (H \\ c) @. A field of /znff = 50 mT 
is also small enough so that the Knight shift [12j and Zee- 
man current [3] effects can be neglected in the analysis 
of the spectra. In the normal state about 6 • 10 6 and in 
the SC state 10 — 20 • 10 6 positrons events were recorded 
for a /iSR time spectrum. 

Typical TF-/iSR time spectra in the normal and SC 
state are displayed in Fig. [1] The temporal oscillations 
and damping of the ^SR asymmetry reflect directly the 
local magnetic field distribution at the muon stopping 
sites. The /j,SR time spectrum consists of two contribu- 
tions: a background signal (Bg) arising from the muons 
stopping in the silver sample holder for the LTF spec- 
trometer or the pressure cell for the GPD spectrome- 
ter and a si gna l arising from the muons stopping in the 
sample (S) [llj. These contributions are clearly seen in 
Fig. Q] At short times the sample contribution domi- 
nates: in the SC state the damping of the signal is en- 
hanced due to the field broadening generated by the vor- 
tex lattice (VL), and the oscillating frequency is reduced 
due to diamagnetic screening. In contrary, for t > 7 fis 
in the normal state and t > 3.5 /xs in the SC state, only 
the signal of the muons stopping in the silver background 
persists. The /iSR time spectra are well described with 
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FIG. 2. (Color online) Local magnetic field distributions 
Ps(B e , T) of CeCoIns at fi H « 50 mT for H \\ c in the nor- 
mal (black squares) and SC (red circles) state obtained with 
the GPD (bottom; p = 0.2 GPa, T ~ 270 mK) and the LTF 
(top; p = GPa, T ~ 19 mK) spectrometers. The blue lines 
are FT of Rs(t). The magenta area represent 0.5Pvl{Bb, T). 
For clarity, the LTF field distributions are shifted vertically. 
A Gaussian apodization (A(t) exp [— 1/2(£/t) 2 ]) of r —A /is 
(5 fjs) was used for the GPD (LTF) spectra. In the inset the 
FT of A(t) (black lines) is shown. The orange dotted lines 
are the FT of R.B g (t). 

the following equation: 

A(t) = A [(l - f s (9))R Bg (t) + fs(9)Rs(t)} (2) 

Here fs(9) denotes the fraction of muons stopping in the 
sample, Aq the initial asymmetry of the signal and Rs(t) 
[Rsgit)] is the sample [background] muon depolarization 
function. fs(9) was determined to be ~ 82% for the LTF 
spectrometer and typically ~ 45% for the GPD spec- 
trometer. Here fs{9) varies each time the pressure cell 
is manipulated (change of p or 9) as the sample position 
relative to the muon beam is modified. In various config- 
urations we recorded a /iSR spectrum after a small field 
increase of 4 mT at low temperatures. Due to pinning the 
shift of field in the sample is less, allowing to determine 
precisely the fraction of muons stopping in the sample. 
The background depolarization function is described by 
a Gaussian field distribution [lT| : 

R Bg (t) = cos(^{B Bg (9,T))t+McM~iyB g (0,T)t 2 /2) 

(3) 

The average background magnetic field (BB g {9,T)) ~ 
50 mT and the standard deviation of the Gaussian field 
distribution aB g {9, T) vary in the SC state since the dia- 
magnetic sample induces a field inhomogeneity in its sur- 
rounding. The initial phase 4>q is constant. 

The sample depolarization function may be written as: 

R s (t)= exp(- 7 24(0,T)i 2 /2) 

x / P VL (B e ,T) cosij^Bgt + <fo) dB e 1 ' 

The presence of a VL gives rise to a local magnetic 
field distribution along the direction of the applied field 
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FIG. 3. (Color online) Pressure dependence of \ a and 7 = 
Ac/Aa in CeCoIns. Red circles are obtained from angular 
dependent spectra, black squares from the temperature de- 
pendent ones. Dashed lines are guides for the eyes. 



Pvl(B s ,T), reflected by the integral in Eq. g]). For 
an extreme type-II superconductor in the London limit 
Pvl(Bo,T) is uniquely determined by an effective pen- 
etration depth A e/ /(0,T) 0[lj|. For the two principal 
magnetic field orientations one has: X e ff(9 — 0°,T) = 
A Q (T) and A e// (0 = 90°, T) = \/ X a {T)X c (T). The 
first factor in Eq. (Q} describes the muon depolariza- 
tion due to additional contributions (cr|(0, T) = cr%{0) + 
, T)) 16]: (i) the nuclear moments [on{9) ~ 



0.5 mT] and (ii) the disorder of the VL [a dVL (8,T) = 
<pyLo(6)\ 2 eff (6,0)/\ 2 ef f{6,T) with a dVL0 (9) « 0.5 mT 
jl7|]. The local magnetic field distribution in the sample 
Ps(Bg, T) can be obtained from the cosine Fourier trans- 
formation (FT) of the experimentally measured A(t), af- 
ter subtraction of RBg{t) (Fig. [2]). 

The angular dependent spectra were analyzed globally 
[Eqs. ([2]) to fli] with the constraint that A e// (6>, T) = 

A a (T) A /cos 2 (0)+ 7 (T) sin 2 (0) with y(T) = X c (T)/X a (T) 
19]. X e ff(8, T — 270 mK) was used to determine the ex- 
act orientation of the sample in the pressure cell (position 
of 9 = 0°). The pressure dependence of the obtained pa- 
rameters A and 7 are shown in Fig. [3] (red circles). The 
analysis of the temperature dependent spectra treated 
globally [HI is also presented (black squares). The two 
data sets give similar results, although different assump- 
tions were made (T or 6 dependence of some parameters 
fixed), demonstrating the reliability of the model. 

The analysis yields at ambient pressure a value of 
X a {T -» OK, fj, H = 50mT) = 350(12) nm. Since in 
the London model the contribution of the vortex core 
is neglected, this value is overestimated [20j]. Taking 
Mo-ff C 2 Orb. — 7.5 T for the orbital upper critical field 
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this correction is only ~ 4% for an applied field of 
fioH = 50 mT and therefore was neglected. In compari- 
son, the first /zSR experiment reported A a (0) « 550 nm 



22]. This experiment was performed in a large mag- 
netic field (J-qH — 0.3 T, and in the analysis an ad- 
ditional field broadening was neglected [16(. This is 
very likely the main reason for the larger value of A a (0). 
Neutron diffraction experiments reported A a (0) between 



247(10) nm [23J and ~ 465 nm |24j. The first value, 



measured in a magnetic field of 2 T, was underestimated 
because Zeeman currents that produce an additional con- 
tribution to the field broadening [25| were neglected. The 
second value was deduced from measurements at 0.5 T. 
Including the correction for the vortex cores [2(J, we 
obtain A a (0) ~ 360 nm in agreement with the present 
value. Surface impedance techniques provided smaller 
A a (0) « 260 nm ;26] and A Q (0) = 281(14) nm These 
experiments were performed in an extremely low mag- 
netic field (< 10/iT) in the Meissner state. 
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FIG. 4. (Color online) Af(0)/Af(T) for the two main orien- 
tations i = a,c for CeCoIns. Interestingly, A^(0)/Af(T) is 
pressure independent in the range investigated. Full sym- 
bol were taken with the GPD spectrometer, empty squares 
with the LTF spectrometer. The empty green circles are ob- 
tained from tunnel-diode oscillator experiments adapted from 
Ref. [2^|. Black lines are power law fits to the data as de- 
scribed in the text. 

The temperature dependence of Af (0)/A^(T) (i = a, c) 
is displayed in Fig. 21 together with a fit of the form 
1 — (T/T c ) n . The pressure evolution of T c was deter- 
mined independently by SQUID magnetometry (Fig.[S]). 
The exponent n was found to be n = 2.17(6) for i = a in 
agreement with Ref. (28|, while n = 1.35(19) for i = c. 
Similar temperature dependences can be obtained us- 
ing AA; = Xi(T) — Xj (0) measured by tunnel-diode os- 
cillator experiments [26j taking A a (0) ~ 336 nm and 
A c (0) ~ 421 nm from this work (green empty circles in 
Fig.|4]). Within precision both A?(0)/Af(T) (i = a, c) are 
pressure independent, suggesting that the gap symmetry 
is unchanged, and the value of the gap to T c ratio is con- 
stant, in agreement with the variation of less than 10% 
obtained by NQR [2^] in this pressure range. 

The pressure independent value of 7 = \J m* / m* ~ 1.3 
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FIG. 5. (Color online) Pressure dependence of 1/A 2 (270 mK) 
(left scale, same symbol convention as in Fig. [3| measured 
by piSR and T c (right scale) determined by SQUID magne- 
tometry of CeCoIns. Dashed lines are parabolic pressure de- 
pendences. Inset: magnetization curves recorded at selected 
pressures and used to determine T c (T c criterion is the onset 
of the transition). 
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FIG. 6. (Color online) Plot of T c versus 1/A 2 (270 mK) 
for CeCoIns under pressure (same symbol convention as in 
Fig. [3]). A linear relation (dashed line T c (X a ) = T c0 + A/X 2 a ) 
with T c0 = 1.88(3) K and A — 4.6(3) ■ 10" 2 K^m 2 is found. 
Inset: field dependence of the magnetic field standard devia- 
tion an (see text for details). 
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(Fig. [3]) is in good agreement with the constant value de- 
duced from the initial slope of the upper critical field un- 
ssure H' c2 {H J_ c,p)/H' c2 {H || c,p) cx m*/m* a « 2 
3fJ. The same experiments indicate that the vari- 
ation of m* with pressure is less than 10% between 0- 
1 GPa with a maximum around 0.5 GPa. Such a small 
variation cannot explain the pressure dependence of 1/ A 2 
plotted in Fig.[S] Therefore, we conclude that, within the 
London model [Eq. (JTJ], ns increases with pressure. Us- 
ing an average value m* ~ 50mrj [5J one obtains from Eq. 
([1]) a change of ns from ns ~ 1.8 to 3.4 carriers per unit 
cell between and 1 GPa. 

In the following we discuss two possible scenarios for 
this strong increase of ns- The first one relies on the 
proposed multigap SC state of CeCoIn 5 Q- The ob- 
served increase of ns with pressure would result from an 
increase of the small gap at 50 mT. Indeed, at ambient 
pressure for fioH « 50 mT the smaller gap is already 
closed To check whether the small gap could open 
under pressure, we probed the field dependence of the to- 
tal magnetic field standard deviation in the sample an- 
Here a% = of + (0.06092<J> /A 2 ) 2 Q is the quadratic 
sum of as previously defined and the magnetic field stan- 
dard deviation generated by the VL ($ is the flux quan- 
tum). For comparison, for the two-gap superconductor 
PrOs4Sbi2 31] different slopes dao/dH are observed be- 
tween the low magnetic field regime with two opened SC 
gaps and the high mag netic field regime where only one 
SC gap is present [32j. In CeCoIns, the fact that the 
magnetic field dependence of an is the same at GPa 
and 0.6 GPa (inset Fig. [6]), strongly suggests that a single 
gap is probed in the full pressure range. 

Another scenario is based on an increase of the Ce va- 
lence (orbital occupancy ny 0.9 [33[ at ambient pres- 
sure). Such a scenario was proposed for the parent com- 



pound CeRhIn 5 where a similar decrease o f A ? is observed 
between p = 2.07 GPa and p = 2.26 GPa [34]. Note that 
valence fluctuations are often associated with SC in Ce 
based heavy fermions (35| . 

An interesting observation is the linear relation 
T c (X a ) = T cQ + A/ XI (270 mK) shown in Fig.© This re- 
lation has some analogy with the Uemura plot [T c (X a ) = 
^uAa(O)] [111 found for underdopped cuprate supercon- 
ductors and other electronically doped unconventional 
superconductors. However, substantial differences exist: 
(i) no proportionality (T c o ^ 0) and (ii) A is about 90 
times smaller than Atj- In addition, pressure affects 1/A 2 . 
much more in CeCoIns than in cuprates [Tlj |. 

In conclusion, we show by TF-/uSR that in CeCoIns 
the magnetic penetration depth (A a ) decreases under 
pressure, while the anisotropy (7 = X c /X a ) and the 
temperature dependence of the penetration depth ratios 
^a(0)/A 2 (T) and A 2 (0)/A 2 (T) are almost unaffected. In 
the range of pressure investigated, a linear dependence 
between T c and 1/A 2 was found. Within the London 
model, the decrease of A Q (270 mK) under pressure cor- 
responds to a doubling of the number density (ns) be- 
tween and 1 GPa, possibly related with the presence of 
a quantum critical point. 
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